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Abstract 

We explore the deeply virtual Compton scattering process at NLO with the emphasis on the 
measurement of the tensor gluon skewed distributions. We recalculate the one-loop coefficient 
function and give two-loop results for anomalous dimensions and exclusive evolution kernels re- 
quired for studying the logarithmic scaling violation. We evaluate weighted cross sections which 
give direct access to these functions on experiment. 
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1 Introduction. 



The studies of the gluonic content of hadrons are comphcated by the fact that there exists no probes 
which couples directly to it. The only way for electroweak bosons to fuse with coloured gluons 
is via quark loops. Obviously this is an effect suppressed in the QCD coupling constant on the 
background of dominating tree level quark contributions. Gluon distributions mix logarithmically 
with quark ones and although they can affect the shape of the latter in a significant way it is still 
extremely complicated to disentangle them in perturbative evolution. An opportunity to measure 
gluons directly could appear provided there is a selection rule which forbids quarks to enter as a 
leading effect in the cross section. 

At leading twist-two level there are three non-local light-ray operators bilinear in gluon fields. 
E.g. in the light-cone gauge where only transversely polarized vector fields do propagate, we have 
the following projections for the Lorentz structure of the composite operator B ^{xi) B^{x2) g^^j^i g^i,i 

where g^^ = g^^ - n^n* - nln^, e^^ = e^upan*n^ and r^^.^^ = i [g^pgt + 9^a9up - 9tu9pa) ■ Here 

and n* are two light-cone vectors such that = n*^ = and nn* = 1. Eq. (|1|) corresponds to 

the Clebsch- Gordon decomposition of the direct product of the two vector representations of the 

Lorentz groupQ (i, |) (g) (i, = (0, 0) © ((1, 0) © (0, 1)) © (1, 1). The first two tensors stand for 

vector and axial sectors and they contribute to the conventional deep inelastic scattering (DIS) 

process on spin-| hadrons. The last one cannot appear there since it requires the photon helicity 

to be flipped by two units. However, it definitely can appear in scattering on spin J > 1 targets as 

was emphasized in Ref. Since these operators belong to the spin-2 representation (1, 1) of the 

Lorentz group, they cannot mix with quark operators. Therefore they can serve as a clean probe of 

the gluonic content in hadrons not contaminated by other effects. Recently this issue was discussed 

in the context of current fragmentation in DIS 0] and in the production of two pions in 77 fusion 

0. Nevertheless this operator can appear in the Compton scattering on the nucleon provided 

it is sandwiched between states with different momenta as was pointed out |0] and elaborated 

in Ref. The corresponding kinematics of the process is known as deeply virtual Compton 

scattering (DVCS) [|, 0, ||. In this note we address the issue of the tensor gluon skewed parton 

distributions (SPD) in DVCS in great detail. We recalculate the one-loop coefficient function, and 

present two-loop anomalous dimensions and exclusive evolution kernels for this sector. We give 

an explicit prediction for weighted cross sections which can be used to extract directly the tensor 

gluon SPD from experimental data. 

^Recall that the representations of the Lorentz group l\_ = 50(3, 1) = 50(4, C)ir « (52.(2, C) (g) 5i(2, C))^^ 

-2 

are labeled by a pair (ji,j2) which are eigenvalues + 1) of the SL{2) Casimir operators J^. 
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Figure 1: One- loop diagrams for the coefficient function. 

2 Leading twist amplitude. 

The hadronic part of the deeply virtual Compton scattering amplitude is defined by the off-forward 
matrix element of the correlator of two electro-magnetic currents sandwiched between states with 
unequal momenta 

T^,{q,P^,P2) = I J dxe'^-SP2\T3^{xl2)3,{-x/2)\P^), (2) 

where q = {qi + '?2)/2 (and the index /i refers to the outgoing real photon with momentum ^2) 
and Pi (P2) is the momentum of incoming (outgoing) nucleon. The leading contribution of the 
light-ray tensor gluon operatoif] 

^(^lA^^i^ 1^2) = G+p{K2n)T^^.p^<l>[K2n, Kin]G^+{Kin), (3) 

whose off-forwards matrix element is parametrized via two SPDs 

Glit,v,A^)^AP-'l^e-'^-{P2fO^^An,-Km 

+E^{t,v,AY-^A^U{P2) [^-Vl,^ U{Pi), (4) 

into the operator product expansion of currents (0) appear at one-loop order (see Fig. 0). The 
traceless symmetric projector r-*- in Eq. (H) possesses the properties T^;po-T^;p/o-' = Tpa-p'a'^ '^^u^pa- = 
Tp^.p^, Tpp-pa- = 0, T^uifiu = 2. The kinematical variables used here and below are introduced as 
uj = = —P ■ q/q^ (generalized Bjorken variable), r] = A ■ q/ P ■ q (skewedness) , in terms of the 
vectors P = Pi + P2, q = |(gi + ^2) and A = P2 - Pi = qi - q2- 

A simple calculation of one- loop diagrams (see Fig. |I]) gives us the following result, with 
restriction of the reality of the final photon being relaxed, i.e. 7^ 0; 

= ~T,Y.Qi Ijt (M, A>(i.,) {1 + . (5) 

^The path-ordered link factor $[2:2,2:1] ensures gauge invariance. 
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which agrees with Ref. |p. The function cr(t, t]) appears from the conversion of gluon field (taken 
in the hght-cone gauge) into strength tensor, = d^B^, and reads, for the fixing of the residual 
gauge symmetry consistent with canonical hamiltonian formalism 0, 

"^''^^ = (t-, + .oKt + ^-^o) - 

Apart from (suppressed) momentum dependence of the coupling constant an additional source 
of scaling violation results from the renormalization of the composite operator (^ which will be 
discussed in the next two sections at one- and two-loop order. 



3 Evolution: anomalous dimensions. 



At leading order in the coupling constant the light-ray operator (|^) obeys the light-cone position 
evolution equation of the form (here and below y = 1 — y) 



d 



din /i^ 
with the kernel 



pO^(«:i, K^)] = ^ I 'dz I ' dy'^'^JC^iy, z)[W{yK, + yn^, zk, + zk^)] (7) 



GG)rT 



[y^z) = CA{y-2 + 



5{z) +Ca{z-2 + 



where /3o = ^TpNf — —Ca is the first coefficient of the QCD /^-function. Fourier transformation of 
this result gives the exclusive evolution kernel known before [[10| , ^ |ll|. The tree level conformal 
invariance of the QCD Lagrangian allows to diagonalize this equation in the basis spanned by 
Gegenbauer polynomials^ 

^Ol = G^,{zd^y-Xu;,.C^Z (5+ (9) 

which form an infinite dimensional representation of the conformal group in the space of bilinear 
operators. Therefore, 

fOj,] = -l-±^^^J,fOj,l with ^Sr = 4C^(^(j + 2)-^(l))+/3o, (10) 



din ij? 



k=l 



the first term in the expansion ^^'jj^ = ^^-fj^^^6jk + {^f ^^ijk^^ + 

In the momentum fraction space the evolution of the SPD can be done making use of orthogonal 
polynomial reconstruction (in the following Gegenbauer polynomials, Cj^'^) of the function from 



its conformal moments according to Ref. |]T3 



k=l 



/ dtClL\{t/r^)G{t,v,Ql), 
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^We drop in what follows the Lorentz indices on the operators. Here d — d+d and V—V — V- 
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Figure 2: Input SPD in (a) and its evolved shape at 4 GeV^ in (b) and at 100 GeV^ in (c) for 

iVmax = 100. 

where formally A/'max = oo. Here Cj^_^i{t) = f^^^j^ (1 — t^)"^ Cj^_^i{t) are adjoint polynomials and 
the re-expansion coefficients Cjk{ri) = (C'f^^i(^) |C'j-i(^i)) are expressed in terms of hypergeometric 
function 2-^1 (?7^). The evolution is demonstrated in Fig. ^ where we have taken an ?7-independent 
inputg, G{t, r^) = 1(1 - 1^), at very low = o.2 GeV^ and evolved it up to = 4 GeV^ (b) and 
g2 = 100 GeV^ (c). 

As a side remark on the solution of the leading order equation for "transversity" let us note 
that even in the case when any arbitrary number of gluons with the same helicity, whose pair- 
wise interaction is described by the kernel (P), are exchanged in the t-channel, one can still 
diagonalize the (LO) multi-particle kernel since the problem admits a large enough number of 
conservation laws to be completely integrable. This can be found by noticing the equivalence of 
the anomalous dimensions ([10|) (which depend on the eigenvalues of SL{2) Casimir operator J ) 



to the Hamiltonian of the exactly solvable one-dimensional XXXs=-3/2 spin chain model ||ri . 

Next we compute the two-loop anomalous dimensions for the tensor gluon operator. To this 
end we use our machinery developed in Ref. [O]. Let us give a brief outline of the method. The 



structure of anomalous dimensions of conformal operators at two-loop (and higher) order reads 

GG T(l) _ GG T(1),D . GGT(1),ND 

where ^''^'yj^^^'^ are the next-to-leading order (NLO) forward anomalous dimension^ |]T^ 

- cl [s,, . 1) (f - .s, {l±l)) - {l±l) . ssu ^ 1) - - f } 



^We factored out the dependence into the gluonic form factor Fg(A^) which for phenomenological estimations 
can be taken equal kt (l — A^/Af^) ^ with Ma = 2.7 GeV and unknown parameter kt which defines the magnitude 
of proton matrix element of the tensor gluonic operator. 

^We thank W. Vogelsang for providing us his result for local anomalous dimensions. 
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Figure 3: IPI Green function for calculation of the renormalization constant Zj^. 

with 

^Ki) = E i s[ (I) = 2^ E m = E ^-=^sm, 
k=i k=i k=i 

and where aj = |[1 + (—!)•']. The non-diagonal elements of anomalous dimension matrix of the 
conformal operators arise due to one-loop breaking of the conformal symmetry and, since the tree 
level conformal invariance leads to diagonal anomalous dimensions, the one-loop special conformal 
anomaly generates two-loop anomalous dimensions. The use of four-dimensional conformal algebra 
provides a relation between the anomalies of dilatation (read anomalous dimensions in question) 
and special conformal transformations via the commutator [D, /C_] = i/C_ which is applied on the 
Green function of elementary fields with conformal operator insertion. To evaluate the commutator 
the knowledge of scale and special conformal Ward identities, with unraveled pattern of symmetry 
breaking for afore mentioned Green function, is indispensable. A careful analysis reveals the result 
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where '^'^7^''°^ are already known LO anomalous dimensions ([l0|) , djk = —aj-k{2k + 3)/(j — k) ■ 
{j + k + 3) for j > k and gjj^ appears as a counterterm required for renormalization of the product 
of two composite operators: integrated anomaly = / d'^x2x-OA{x) = J d^xx^Z^ {Gfj.u)'^ in the 
trace of energy-momentum tensor and a conformal operator. The structure of the counterterms 
has been established using the form of counterterms for differential vertex operator insertions and 
found to be 



^[04x)]fOJ,] = ^[OA{xfOJ,] - 6^'\x) E {Za}^^ M] - ^^d^S^'^x) E {^a},, " • • • 

k=0 k=0 

where X = ZgZ^^^ is expressed in terms of the renormalization constants of the gluon wave 
function Z3 and the coupling Zg. 

To determine the unknown renormalization matrix Z^ it proves convenient to work in the light- 
cone position formalism and calculate IPI Green function {[O a{x)]['^0'^ {ni, K2)]B ^{xi) B^{x2)) a.mp 
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shown in Fig. |^ (explicit one-loop graphs can be found in Ref. [|T2]). We write the divergent part of 
the operator product K2)] via the following relation in leading order of the coupling 

constant: 



Ztc e Jo Jo I 

+ """"iC^iy^z) Jd''x2x.Bl{x)j^j^^fO''{yK,+yK2,zK^ + zK2)]y{lQ) 

For the evolution kernel "^""^/C^ we find immediately from Eq. (|15]) that '^'^IC^{y, z) = —i^6{y)6{z). 
Explicit calculation of the amputated Green function of conformal anomaly and conformal operator 
gives 

^^IC'Ay, z) = -t2{K, + «:2)^^/C^(l/, z) - z{k, + K,)f3o6{y)6{z) + ^^/C^(y, z), 
where ^^lC^{y, z) is the LO evolution kernel (H) and 



(17) 



'ICl{y,z) = 2 



k 



2+ 

+ 2^ 





T 




.z. 




'2 




y. 



5{z)- 



5{y) + {2y - mz)^ 
5{z) + {2z - my)] . 



The evaluation of the conformal moments of (|1^) [multiplied by 1/2 ■ (j — A;) ■ (j + A; + 3)] gives us 

(3 + 2fc) 



X 1 2Ajk + {Ajk - ^{j + 2) + V^(l)) 



(19) 



r(j + 4)r(fc) 

r(j)r(A; + 4) 

where we have introduced the matrix A with its elements defined by 



A 



■jk 



'j + k + 4^ 



- k^ 



+ 2^(j-A;)-7/;(j + 2)-V^(l), 



(20) 



and 6jk = 1 for j > k and zero otherwise. Note that this result differs from our previous ones for 
chiral even operators only by a pure rational function 

m 



The solution of two-loop evolution equation is straightforward and can be found in Ref. |13 



(21) 



4 Evolution: two-loop kernel. 

In this section we give our results for the two-loop exclusive evolution kernel reconstructed from 
the known non-diagonal anomalous dimensions found above and two-loop splitting function of 
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Ref. |T^. The structure of ER-BL kernel to O (af) accuracy reads 



GGyT 



2tt 



GG 



2^ 



GG 



1GG„T(1) 
2 



where the +-prescription is defined according to 



[Vix, y)] , = V{x, y)-5{x-y) / dzV{z, y). 



7r'^5(x - y) , (22) 



(23) 



Note that this prescription is not in one-to-one correspondence with the +-prescription in the 
forward case. The LO kernel reads 0, ||, 0] 



GG 



V^^'^ = CA9{y-xf^f{x,y) + 



X 

y 



X 

y 



with 



GGrT 



r 



X 



1 



y^y-x 



(24) 



To obtain the NLO correction we use the approach described in [jT5[. Since the tensor gluon sector 
is almost analogous to handle as the vector and axial ones, we only mention the differences to 
these cases. The structure of two-loop anomalous dimensions (0) implies the following form of 
the kernel 



GGyTil) ^ _GGyT 



GGyT(0) + ^ J _ 



_^ggqT^ggjjT^ (25) 



where the commutator stands for [A0 B]{x,y) = Jq dz {A{x, z)B{z,y) — B{x, z)A{z,y)}. The 

off-diagonal part is contained in the first two convolutions on the r.h.s. of this equation and they 
are known exactly. The conformal moments of the dotted kernel, 

„2 



GG 



V^^(°) = CA9{y - x) 



X 



1 , X 

In - 



y y-x y 



X 

y 



X 

y 



(26) 



are proportional to the commutator of the d-matrix with the LO anomalous dimensions. The 
(yf-matrix corresponds to the kernel 



GGgT 



CAe{y-x) 



2^- 



(1 - f ) 

y-x 















[y- 


■^yj 



(27) 



which contains a part (second term) of the chiral even (odd) case. The difference is a rational 

function (first term) which can be restored from the GQ-channel result by applying appropriate 

convolutions in order to trade the denominator + 1) ■ [j + 2) of the GQ conformal moments 

of x/y to go over into the GG sector^, i.e. 1/j ■ (j + 1) ■ (j + 2) ■ (j + 3). 

^This can be achieved by convolution with the so-called c-kernel jl5| in QQ channels and differentiation with 
respect to y, i.e. (Py^y^/dy^. Finally, we removed a pure diagonal piece. 
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Now the problem is reduced to the reconstruction of diagonal parts. The third term reads 
g^gT = J^eiy - x) y) + A^«/.^(x, y)] + e{y - x) y) + A^«/i^(x, y)] ] 



+ 



(28) 



X X \ 

and arises from the crossed-ladder diagram and can be obtained by means of A/" = 1 supersymmetry 
from the one in the quark sector. The particular contributions are 

2x 2x _ 2x 
— ^ In X In 



GGf^T 



2^^/^lnxlni/ - 2^^f ' [Li2(x) + Li2(y)] , A^^/i 



, GGrT 



GGuT 



GGu 



X 



2Li2 ( 1 - - 1 + In" y 



y.y y.y yy. 

+ 2^^/^ [Li2(^) - \nx\ny] + 2 ^^/^Li2(a;) 



(29) 



Here Li2(x) 



/o^ Y ln(l — t) is the Euler dilogarithm and we used the following shorthand 



notation '^'^/^ = '^^f'^[x,y). The remaining diagonal part D has a simple representation in terms 
of LO kernels and can be obtained by taking the forward limit and comparison with the known 
DGLAP kernel [IT3. Consequent restoration of diagonal ER-BL kernels from the splitting functions 
is straightforward and gives 



GGjjT 



13 



GG a , f GG c 

3 



GG^a 



3gg^« 



6 4 



GG^c 



(30) 
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Here the h kernel coincides (with colour factor being dropped) with the LO kernel ( ^If ) and the a 
and c kernels having the structure f*(a;, y) = f^6{y — x) + f^6{x — y) are defined by the functions 
GGja ^ x'^/y'^ and '^^f = x'^{2xy + y - x)/y'^. Eqs. (|25l-|0l) is our final result for the NLO 
corrections to the gluon "transversity" evolution kernel. 



5 Cross sections. 

Now we are in a position to study the cross section for electroproduction of real photon where 
the gluonic SPD, whose perturbative properties we studied in detail in the previous sections, 
can be measured. For DVCS the skewedness r] and generalized Bjorken variable ^ are propor- 
tional Tj = (j- ~^ • '^^^ differential cross section with unpolarized lepton beam and 
unpolarized target, in variables y = Pi ■ qi/Pi ■ k, Q"^ = —qf and x = Q^/(2Pi ■ qi) with 
^ = x(l + ^) (2- x + x^^ , reads 



\ — 1 /2 

da a^xy ( AM'^x\ 



dxdyd\/S?\d<Pr 87rQ2l Q2 



2 



(31) 
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Here (j)r is the azimuthal angle between the lepton and proton scattering planes in the rest frame of 
the target. In the following we will be interested only in the interference term |Xp = ^h^vcs + 
^vcs^H between the DVCS, Tdvcs, and Bethe-Heitler amplitude 7bh since it provides a unique 
opportunity to extract the real/imaginary part of the gluonic amplitudes 



)Nf ^ 
n ^ J-i 



An explicit calculation gives for unpolarized settings 
f±l) 



(32) 



X 2 

«3 



P-q 



where 



ri = 2(Fi + F2)Re Hg + 



4M2 ^ , 



1 



(33) 



(34) 



and +(— ) sign stands for electron (positron) beam. In consequent evaluation of this expression 
we perform an expansion in 1/ and keep only the first non- vanishing contribution. We form 
the charge asymmetry in order to extract the interference term from the total cross section. Since 
the double helicity flip amplitude presently considered is uniquely proportional to cos (30,.) (while 
all other contributions to the non-polarized cross section |X|2 enter with cos(n0,.) and n = 1,2 



I6| one can isolate this purely gluonic contributions from the effects of quark SPD by forming 

) = cos (30r) 



an appropriate weighted cross section / d(f)rw{(j)r)(y{(f)r) Namely, choosing w[ 
we extract tensor gluon SPD from the unpolarized cross section, 

~ 3/2 



27r 



27r Jo 



where dM. 



d(f)r COS {3(j)r 



d'^a — da 



16a 



A2 



7(l-x)3(l-y) 



xy (2 — x) 



1 



A 



A2 



T2dM, 



(35) 



cfixy 



(1 + 



22 



-1/2 



dxdyd\t^\ and A 



2 

min 



M2x2/(l-x + xM2/Q2). Note that 
-A2/ Q}\ while the combination 



the correction to the structure function is suppressed by O 
Ti is down by O (hP/^f^S^) relative to M2r2. 

The imaginary part of the amplitudes and £q can be accessed by means of single spin 
asymmetry. The polarization of the lepton beam does not induce contribution of gluon "transver- 
sity" distributions into the cross section in agreement with However, once the proton beam 
is longitudinally polarized we can access this asymmetry due to the interference term since the 
latter has genuine sin (30^) azimuthal angle dependence. Thus the SPD can be extracted via the 
following weighted cross section 

3/2 

lQ\l-=z^^ ,J ^.."' 11- =^ I fdM, (36) 



1 /•2'^ 
2^ Jo 



■ sin (30,. 



d^ 



d^ 



Q2 xy (2 — a;)2 



A2 
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where in d'^a^ the proton is polarized along the positron beam and 

Let us stress that the results (^) and (|36| ) are valid to accuracy and further expansion 

terms from other structure functions can mimic the cos (30) / sin (30) behaviour and contaminate 
the double helicity flip cross sections. 

6 Conclusion. 

In this paper we have presented weighted real photon electroproduction cross sections which can be 
used as a probe for the magnitude of gluon content in the nucleon. The opportunity to extract the 
real and imaginary parts of these amplitudes is offered by the scattering of the unpolarized lepton 
beam on the unpolarized and longitudinally polarized nucleon targets, respectively. We have given 
as well the formulae for NLO coefficient function as well as two-loop anomalous dimensions of the 
conformal operators and momentum fraction exclusive kernels thus completing the set of results 
required for study of the scaling violation for all twist-two SPDs at NLO. 

This work was supported by DFG and BMBF (D.M.). 
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